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Abstract 

Usually, the Kohlrausch (stretched exponential) function is employed to fit the lu- 
minescence decays. In this work we propose to use the Tsallis distribution as an 
alternative to describe them. We show that the curves of the luminescence decay 
obtained from the Tsallis distribution are close to those ones obtained from the 
stretched exponential. Further, we show that our result can fit well the data of 
porous silicon at low temperature and simulation result of the trapping controlled 
luminescence model. 
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1 Introduction 



The Kohlrausch (stretched exponential) function has been largely employed to 
describe luminescence decays of different materials and in different time scales 
(see p1[2][3]Hll5]l6|7] . and references therein), and it has the following form: 
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with < P < 1. We note that the value of parameter (3 in the range < 
P < 1 makes I{t) decays slower than the exponential function, whereas for 
P > 1 the function ([T]) is called compressed exponential and it has short-tails, 
consequently, it decays faster than the exponential function. 

For the general consideration of the study of luminescence decay I{t), it is 
usually considered the distribution of rate constants H{k) [4J connected with 
I{t) by the Laplace transform given by 



with /(O) = 1. In particular, H{k) must be nonnegative for all > in order 
to be considered as a distribution function. Moreover, the function H{k) is 
normalized, i.e., H{k)dk = 1. We note that H{k) may have a large change 
for a small change in I{t). In this way, the precision of the experimental data 
will be important for the choice of the distribution of rate constants H{k). 

Modification and generalization of the stretched exponential function as a 
decay function has been recently considered, for instance, a decay function 
unifying the modified stretched exponential and Becquerel decay laws [4ji5jt6j. 

In this work we consider a simple distribution of rate constants H{k) based on 
the Tsallis distribution pi9] . Our idea of employing this distribution is due to 
its simplicity and it has been applied to a variety of natural systems (see, for 
instance, p^lllll2p3] ). We shall show that the luminescence decay function 
generated by the Tsallis distribution can be useful to describe experimental 
data. To do so we shall fit the data of porous silicon at low temperature |3J and 
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the simulation result of the trapping controlled luminescence model described 
in [l] with our luminescence decay function. 



2 Distribution of rate constants H{k) given by the Tsallis distribu- 
tion 

For the luminescence decay I{t) given by Eq. ([1]) the distribution H[k) can 
be expressed in an integral representation [4J as follows: 



Hsik) = — duexp 
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We see that this last expression is quite complicated. Now we consider the 
Tsallis distribution for H{k), and it is given by 



Hg{k) = B [1 - {1 - q)ak]^^ (4) 

where a > and < g < 2, and i? is a normalization factor. The parameter a 
has the dimension of time. For q 1 we recover from Eq. (jlj) the exponential 
function which is connected with a particular case of the Becquerel decay 
function ^ given by 

m = -^—T . (5) 

We note that Hq{k) has a cutoff for g < 1, i.e., the term [1 — (1 — q)ak]~ 
is replaced by zero when [1 — (1 — q)ak] < 0, then the normalization factor is 



3 



equal to 

B = a{2- q) (6) 

where < g < 2. The luminescence decay I{t) is obtained by substituting Eq. 
(jl]) into Eq. ([2]), and we arrive at 

2 — a t /"i 1 tu 
^qi^) = -. / C^M It 1-9 6 "(1-9) , < g < 1, (7) 

and 

2 n t 1 tu 

hi.'t) = -e'^^i^ / du wT^e^^ti^ , 1 < g < 2. (8) 

q — I Jl 

In order to see the behavior of Iq{t) we plot some curves of Iq{t) and I pit) 
with typical values of a, q and (3 which are shown in Fig. 1. All the curves 
decay monotonically with time, and two of the curves of Iq{t) have fat tails. 
We note that the lowest curve of Ig{t) in the figure is very close to the curve 
of Ifsit); This means that any experimental data fitted with these curves must 
have a good level of precision in order to choose which one of the curves is 
preferable. 

For application of our result we consider the experimental data of porous 
silicon at low temperature given in [3] . In Fig. 2 we compare the best fit of the 
data by using the stretched exponential and Ig{t). We see that both curves 
can fit the data (which are not shown in our figure) very well. Then, in this 
case, the experimental data do not offer us sufficient precision to discard one 
of the curves. We note that the behavior of Ig(t) is not the same of //^(t), i.e., 
I pit) describes a straight line, whereas Iq{t) does not. 

Another application of our result is to consider the model of trapping con- 
trolled luminescence given in Ref. [1]. The process during the excitation is 
described by the following equations: 
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cLti 

= X - R{M -m)n^ (9) 
— = -A^mn^ + R{M -m)ny (10) 

n 77 

— = An{N-n)n, (11) 
dric dm drii, dn 

^ " IT^ ^ ~ dt ' ^ ' 

where ric and are the instantaneous concentrations of electrons in the con- 
duction band and holes in the valence band, respectively; An is the retrapping 
coefficient, is the recombination coefficient and R is the trapping coeffi- 
cient of free holes during the excitation; and M are the concentrations of 
the traps and recombination centers, whereas n and m are their respective 
instantaneous occupancies. Finally x is the rate of production of electrons 
and holes by the excitation irradiation. The luminescence emission intensity 
is calculated by the rate of electron-hole recombination given by 

-Ammuc . (13) 



dt 



In Fig. 3 shows a replotting of the simulation result [Ij and the best fitted 
stretched exponential function with the parameters x = 10^^m~'^s~^, A^ = 
10-^'^mh-\ An = 10~9m3s-\ R = 10-^Ws-\ = lO^^m-^, M = lO^^m'^. 
We note that the agreement is not very good. In order to enhance the ffiting 
result the decay curve has been separated into two parts [T]: The values have 
been ffited separately for the ffist microsecond and for the period of time from 
t = 2 to 10/is. In Fig. 4 shows the simulation result fitted by our result; we see 
that the agreement is excellent and the decay curve has not been separated. 
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3 Conclusion 



In this work we liave considered tlie Tsallis distribution as a distribution of 
rate constants H{k). From Eq. ([2]) we have obtained the luminescence de- 
cay function Ig{t) for fitting the luminescence decays. We have shown that 
Iq{t) may have fat tails and its behavior looks like the stretched exponential 
function. We have also demonstrated that Iq(t) can be useful to describe ex- 
perimental data. In the case of the trapping controlled luminescence model 
the fitted result by Iq{t) is superior than the best fitted stretched exponential 
function. It can be seen that, Fig. 3, the deviation of the best fitted stretched 
exponential function from the exponential one is not small f3 = 0.45. On the 
other hand, Fig. 4 shows the simulation result fitted by Iq{t) with q = 1.04 
which presents a small deviation of Hq{k) from the exponential function, then 
Hq{k) is close to the Becquerel distribution of rate constants Q;exp(— A;a); 
however, the parameter a = 13.0321/is has a value close to r = 11.8/is. 
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Figure Captions 



Fig. 1 - Plots of Iq{t) and I/sit) given by Eqs. ([T]) and ([HD in arbitrary units. 
The dotted lines correspond to Iq{t), whereas the solid lines correspond to 
Ipit) with typical parameter values: from top to bottom, (3 = 0.6, r = 4.56, 
a = (1.64)2 , g = 1 + (0.28)2 ; p = 0.65, r = 3.07, a = (1.4)^ ,q = l + (0.35)^ 
; P = 0.45, T = 0.68, a = (0.695)^ ,q=l + (0.275)^. 

Fig. 2 - The best fit to Eq. ([1]) (solid line) for porous silicon at low temperature 
described in [3] with /3 = 0.75 and r = 2.4ms. The dotted line is the plot of 
Eq. ([HD with q=l + (0.25)2 and a = (1.235)2ms. 

Fig. 3 - Replotting of the simulation result calculated in [Ij and the best fitted 
stretched exponential function with P = 0.45 and r = 11.8 /is. 

Fig. 4 - Plots of the simulation result calculated in [Ij (solid line) and Iq{t) 
(dotted line) with q=l + (0.2)2 ^^^^ ^ ^ (3.61)2^s. 
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